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The extrapolation ol truncated perturbative series to the 
large coupling regime is not reliable since the accuracy of the 
calculations cannot be improved systematically by a finite set 
of higher order terms. Instead, a nonperturbative approach 
is proposed here to describe the thermodynamics of strongly 
coupled plasmas. Within the (f> 4 theory, it is derived from 
the Luttinger-Ward theorem and amounts to a quasiparticle 
description. The approach is extended to the quark-gluon 
plasma. Even close to the confinement temperature it is found 
to be in a good agreement with lattice data of several gauge 
systems. The model is applied to predict the equation of state 
of the quark-gluon plasma with realistic quark masses. 



I. INTRODUCTION 

Perturbation theory is a well established framework for 
approximative analytical calculations in quantum field 
theory. In numerous examples perturbative calculations 
have been shown to be successful in the description of 
processes and systems characterized by a weak coupling. 
For larger values of the coupling strength, however, ex- 
trapolated perturbative results do not yield appropriate 
approximations. So, e.g., the perturbative series of the 
deconfined QCD thermodynamical potential f2, which is 
known up to the order 0(g 5 ) jy, shows even for moder- 
ate values of the coupling, say g ^ i , a strongly fluctuat- 
ing behavior increasing with the order of the expansion, 
whereas lattice gauge simulations HH exhibit a rather 
smooth dependency fl{g). 

This feature can be understood qualitatively by the 
asymptotic expansions Z n (g 2 ) — J2m Zm (9 2 ) m °f the 
function Z(g 2 ) — f_ dx exp { — i x 2 — g 2 x 4 }. This sim- 
ple example Q mimics the perturbative and full repre- 
sentation, respectively, of the partition function of an 
interacting field theory in zero dimension. (In fact, the 
coefficients z m count the number of diagrams contribut- 
ing to the thermodynamical potential of the scalar </> 4 
theory considered below.) Although Z(g 2 ) defined for 
complex g 2 possesses a cut along the negative axis (thus 
the series expansion is divergent), Z n still contains use- 
ful information as to be expected for perturbative results. 
The error of the approximation can easily be estimated 
to be \Z(g 2 ) — Z n (g 2 ) \ ~ (ng 2 ) n , hence the optimal order 
to truncate the series expansion is n* ~ l/ff 2 - 

These considerations generalized suggest that for the 



(at least qualitative) description of strongly coupled sys- 
tems suitable lower-order perturbative approximations 
might be more appropriate than elaborated high-order 
calculations. For the QCD thermodynamical potential, 
this leading- order conjecture is also supported by the ar- 
gument of minimal renormalization scale dependence of 
the truncated perturbative expansion. With this assump- 
tion for the regime of larger coupling strength, in Section 
the thermodynamical potential of the scalar 4 theory 
is calculated in a selfconsistent nonperturbative approxi- 
mation H . The emerging picture motivates a quasiparti- 
cle description of the thermodynamics of gauge systems 
as discussed in Section [II. The model is compared to 



available lattice data, and a prediction of the equation of 
state of the quark-gluon plasma (QGP) including strange 
flavor is given. 



II. TADPOLE RESUMMED 4 THEORY 

It is well known that the (naive) perturbative calcula- 
tion of the thermodynamical potential f2 of the massless 
scalar theory with the bare Lagrangian 



C 



(i) 



is plagued with infrared divergencies §. These ther- 
mal divergencies can formally be treated by suitably re- 
arranging the Lagrangian, e. g. by adding and subtract- 
ing a mass term to C, yielding an equivalent theory with 
free massive propagators and an additional 2-point in- 
teraction. With such a reorganization of the perturba- 
tion theory, the expansion of £1 in the coupling strength, 
then well-defined up to 5th order, has been calculated 
0. However, this series expansion exhibits a similar bad 
convergence as that of the QCD thermodynamical po- 
tential when being extrapolated to larger values of the 
coupling. As an alternative way to the perturbative ex- 
pansion around the (modified) interaction-free limit, the 
Luttinger-Ward approach, originally derived for nonrela- 
tivistic fermion systems || but readily generalized ||, is 
proposed here. Formulated in terms of full propagators, 
this formalism appears to be a more adequate description 
than perturbative calculations using free Green's func- 
tions. In the tadpole approximation considered below 
it leads to propagators with a mass-like selfenergy pro- 
viding thus a less formal motivation of the perturbative 
reorganization technique. 
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The Luttinger-Ward formulation of the thermodynam- 
ical potential of the scalar theory reads 



Of 



where A = (Aq" 1 — II) _1 denotes the exact propagator 
and II is the exact selfenergy which can be decomposed 
into skeleton (i. e. 2-particle irreducible) contributions 
11^ of order m. The sum-integral notation is specified 
below. The expression (|^) obeys the fundamental condi- 
tion g(J 
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which functionally relates microscopic properties (II) to 
macroscopic quantities (f2) ensuring thermodynamical 
selfconsistency. Moreover, the representation (||) is an 
expedient starting point for a systematic (also referred 
to as symmetry conserving) approximation scheme. Ap- 
proximating fl' by the order-n truncated series Q,' n and 
the selfenergy by II n = J2m yields an approximative 
expression 51„ for the potential which strictly fulfills (||) 
& 

Relying in the following on the conjecture stated in 
Section M, f2' is approximated in leading-loop order by 



where double lines represent the full 



propagator Ai = (A 1 — TIi) 1 within the approximation 
considered. The corresponding selfenergy is given by 



IIi = 12 
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This implicit equation for IT amounts to the summation 
of all super daisy selfenergy diagrams Q of the (naive) 
perturbation theory. Due to the tadpole interaction pro- 
cess, IIi is real and momentum independent. Thus, Ai 
takes the form A m 2(P) = (P 2 — to 2 ) -1 of a free propa- 
gator with a mass term so the sum-integration in (0) is 
elementary, though divergent. It can be regularized in 
4 — 2e dimensions using in the imaginary time formalism 



T 



47T 



T 



E 

Po 



2 c 



P 



(27T) 



3-2e 



, po = i 27rTn, with 



7 being Eulcr's number and the MS renormalization 
scale, 
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(e). (5) 



Here the notation lj p = (to 2 + p 2 ) 1 / 2 and the Bose func- 
tion 71b{x) = [exp(x) — are used. The first expres- 
sion of the right hand side of equation (H) corresponds 



to the vacuum selfenergy term of a massive scalar par- 
ticle. It contains a divergent contribution which in the 
present case cannot be compensated by a counter term 
since it is temperature dependent. However, both the 
divergent and the scale dependent term are absorbed 
when expressing go by the physical coupling g in equa- 
tion (Q). To be consistent with the tadpole topology 
class for the selfenergy, g is defined as the truncated vac- 



uum scattering amplitude 



at the momentum scale s: g 2 = g$ — ^ <? 2 <?o L(s) with 
L(s) = {An)- 2 [e- 1 - ln(s//2 2 ) + 2 + 0(e)] . With the 
identification Hi = m 2 and choosing </s = eT, the renor- 
malized equation (Q) reads at e — * 
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The solution m 2 (g) of this gap equation, when expanded 
in <?, coincides with the perturbative selfenergy up to 
0(g 3 ). For larger values of g, in contrast to the pertur- 
bative result in the various known orders, it has a smooth 
behavior. 

In the tadpole approximation, the thermodynamic po- 
tential (0) is given by 



ln(-Ar 1 ) + iTO 2 A 1 
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where @ and d m i ln(— Aj[ x ) = -Ai were used for evalu- 
ating the sum-integral over ln(— A] -1 ) up to an irrelevant 
integration constant. As to be expected for a consistent 
approximation, the divergent and renormalization scale 
dependent terms from the sum-integrals exactly cancel 
each other in In analogy to the approximation Hi of 
the selfenergy, sums up all perturbative super daisy 
diagrams of the thermodynamical potential. 

It is emphasized that in equation (|7|) the first contribu- 
tion is the thermodynamical potential of an ideal gas of 
quasiparticles with mass m(T) as determined by the gap 
equation (^). The remaining interaction contribution, 
when calculating the entropy S± = —dQ\/dT, compen- 
sates the temperature derivative terms of the quasipar- 
ticle mass. Thus, according to (|J), Si is given by the 
entropy of the ideal massive gas. 

Figure |l| shows the pressure p\ = — Sli/V^ as a function 
of the coupling strength. In the range of the coupling 
considered, the positive interaction pressure turns out to 
be a correction to the ideal gas contribution pid {T, to) = 
T/(2tt 2 ) / dpp 2 \n(l - cxp{-w p /T}) at the 10% level. 
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FIG. 1. The pressure p\ and the ideal gas contribu- 
tion J5id as functions of the coupling strength g in units 
of Psb = 7r 2 /90T 4 . For comparison, the leading and 
next-to-leading order perturbative results are shown as well. 

For large values of g the pressure p\ does not deviate 
as much from the free limit as the naive extrapolation of 
perturbative results might indicate. This general feature 
is confirmed by other nonperturbative approaches. So the 
approximation f2i agrees with the result ]T^] derived by 
considering the 0(N) scalar theory in the limit N — > oo. 
Without an expansion in g, the reorganized (screened) 
perturbation theory fll4| , where a gap equation has to 
be taken as an 'external information', leads to similar 
results. 

Compared to these approaches, the Luttinger-Ward 
ansatz is favored here because, by construction, it yields 
symmetry conserving approximations. Furthermore, ex- 
pressing the thermodynamical potential in terms of full 
propagators, the formulation (^|) presumably guarantees 
a well-organized asymptotic expansion ]l5[ |. Finally, it 
is at least conceptionally straightforward to extend the 
selfconsistent resummation scheme to the next-to lead- 
ing order where the functional fl' is approximated by 



n f 2 = n[ 





Due to the rising sun contribution 



the selfenergy then becomes momentum de- 
pendent and complex thus the quasiparticles acquire a 
finite spectral width, a fact which requires more sophis- 
ticated techniques as those utilized here. 



III. A QUASIPARTICLE MODEL FOR THE QGP 

In QCD the difficulties mentioned at the end of the 
last section arise already at the 1-loop level. There- 
fore, in the following a more phenomenological approach, 
though closely related to the Luttinger-Ward formalism, 
is applied to describe the strongly coupled QGP. Temper- 
atures closely above the confinement temperature T c ~ 
200 MeV, as they are possibly attainable in heavy-ion col- 
lisions, thereby require the appropriate treatment of the 



strange flavor. With a mass m,Q S ~ 150 MeV the strange 
quarks represent a relevant non-light degree of freedom 
while the heavier flavors are suppressed. 

The plasma is characterized by broken Lorentz sym- 
metry which leads to additional partonic excitations not 
existent in the vacuum. The effect of the symmetry 
breaking is yet marginal for excitations with momen- 
tum p ^ T. Therefore, the hard 'medium-bound' (lon- 
gitudinal gluonic and helicity-flipped fermionic) modes 
are essentially not populated while the persisting exci- 
tations bear resemblance to the physical vacuum states. 
These thermodynamically relevant excitation are approx- 
imately described by a dispersion law ui 2 (p) = m 2 + p 2 



with the effective gluon and quark masses given by [p|]lG|] 



m 9 = \h [N c +-rH )g 2 T 2 



N 2 - 1 



m 0q + \l m 2 q + ^ g 2 T 2 



(8) 



These gauge invariant masses are generated dynamically 
and thus not in conflict with two propagating gluon ex- 
citations and with chiral symmetry for vanishing current 
quark masses ma q . Within the mass parametrization 
(g) of interaction, sub-leading effects in g like damping 
are neglected so the approach outlined below amounts 
to a partial resummation of relevant contributions simi- 
lar to the Luttinger-Ward formalism. The expansion in 
m/T ~ g yields the leading order perturbative QCD re- 
sults 0. 

The pressure of this system of effective quasiparticles 
with temperature dependent masses TOj(T) is given rem- 
iniscently to @ by || 



PqpCn = I>d(r, m t (T)) - B{T) . 



(9) 



The function B(T) resembles the interaction term in (0). 
It is not an independent quantity. From the thermo- 
dynamical relation (|^), which in the present framework 
takes the form dp qp /dmi = (implying for chiral quarks 



0, thus the model respects chiral symmetry 



restoration in the plasma), B is determined by 
dB dPid drrii 



dm,; dT 



(10) 



As a special case, constant effective masses lead to 
B(T) = const and result in a bag model equation of 
state. 

The entropy density s = dp/dT tests the occupation of 
the phase space, consequently s qp (T) = Sid(2\m$) is 
insensitive to B(T). This function has the obvious inter- 
pretation of the ground state energy of the quasiparticle 
system since the energy density e = sT — p reads here 
e qp (T)=J2 i e id (T,m l (T)) + B(T). 
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To account for nonperturbative behavior near T c , g 2 
in (0) is specified by the effective coupling 



G 2 (T) 



48tt 2 



(11) 



(f f N c - 2 n { ) ln[(T + T S )/(T C /A)] 2 

where the parameter T s /T c is a phenomenological regu- 
lator. For large temperatures, G 2 approaches the 1-loop 
running coupling with T c /A ~ Aqcd- 

Equations (g) - (0) constitute the effective quasipar- 
ticle model which selfconsistently 'continues' the lead- 
ing order perturbative QCD thermodynamics. Corrobo- 
rating the leading-order conjecture, it quantitatively re- 
produces thermodynamical lattice gauge data |H,p| by 
suitably adjusting the parameters A and T s /T c ]17 ,|l9|| - 
Hereby it turns out that fixing the integration constant 
in ( |l0| ) at the confinement temperature yields positive 
values B a = J2iPid(T c ,mi) - pi a tticc(T c ). This fact sug- 
gests a bag constant interpretation of Bo. In both the 
pure gauge system and the flavored plasma, however, at 
a temperature T ~ 2T C the 'bag' function B(T) becomes 
negative due to the asymptotic behavior of the coupling 
which is opposed to the scalar model. 

Motivated by the comparison to available lattice data, 
the quasiparticle model is now applied to predict the 
equation of state of the QGP with realistic quark masses 
which so far cannot be simulated numerically. The con- 
finement temperature and the bag constant are chosen 
to be T c = 170 MeV and B = (180 MeV) 4 . Equat- 
ing the pressure of the confined phase (described by a 
hadron resonance gas) and p qp at T c to fulfill Gibbs' con- 
dition relates the remaining two parameters A and T s /T c . 
In Figure | the predicted QGP pressure, energy density 
and the function B(T) are shown for the large variation 
3 < A < 10 of the free parameter. 




TIT C 

FIG. 2. The equation of state of the QGP with massive 
strange flavor as predicted by the quasiparticle model, and 
the function B(T). Decreasing values of the model parameter 
3 < A < 10 are indicated by the arrow; the dotted short 
line marks the asymptotically free limit. The dashed lines 
represent pressure and energy density of the confined QCD 
phase approximated by a hadron resonance gas. 



rather insensitive to the choice of A, T c and Bq- Re- 
markably, for temperatures T ^ 2T C the energy density 
shows a saturation-like behavior at some 80-90% of the 
free limit. It is worthwhile pointing out that the strange 
flavor yields an important contribution to thermodynam- 
ical quantities because the effective strange quark mass 
is yet smaller than the thermal gluon mass. 

In summary, a QGP quasiparticle model motivated by 
the Luttinger-Ward approach to 4> A thermodynamics is 
presented. Being successfully tested on lattice data of 
various gauge systems, it provides a realistic prediction 
of the equation of state of the QGP with strangeness. 
This allows further quantitative estimates of relevant ob- 
servables in relativistic heavy-ion collisions. 
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